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Abstract 

In ultrarelativistic heavy-ion collisions at y/s > 20 (120) A-GeV a copious production 
of charm (bottom) production sets in which, via correlated semileptonic DD (BB) de- 
cays, gives rise to a dilepton yield at invariant mass M ~ 2 - 3 GeV in excess of the 
Drell-Yan yield and the thermal dilepton signal from deconfined matter as well. We show 
that appropriate single-electron transverse momentum cuts (suitable for ALICE at LHC) 
cause a threshold like behavior of the dilepton spectra from heavy-quark meson decays 
and the Drell-Yan process and can allow to observe a thermal dilepton signal from hot 
deconfined matter. 
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1 Introduction 



The current relativistic heavy-ion experiments at CERN-SPS and planned future exper- 
iments at BNL-RHIC and CERN-LHC are aimed at a study of the properties of highly 
excited, deconfmed matter. The production of dileptons with intermediate invariant mass 
in the so-called continuum region between and J / if) are widely considered as one of the 
main tools to measure in such experiments the "initial" temperature and other ther- 
modynamical characteristics of the produced matter [[I], Dilepton measurements are 
envisaged in particular with the PHENIX and ALICE detector facilities at RHIC and 
LHC, respectively. The thermal dilepton signal in the intermediate invariant mass region 
faces a serious background problem connected with hard initial quark - anti-quark anni- 
hilation into dileptons (i.e. the Drell-Yan production). The general expectation is that 
with increasing beam energy the maximum temperature of matter rises and consequently 
the thermal dilepton yield grows faster than the Drell-Yan yield. However with increasing 
beam energy also a copious production of heavy quarks, resulting in hard initial collisions 
of partons too, sets in. As a consequence, the correlated semileptonic decays of open 
charm and bottom mesons yield a dilepton rate exceeding the one from the Drell-Yan 
process. This can be understood by the following simple arguments. Drell-Yan pairs are 
produced dominantly via qq annihilation, while the open charm and bottom production 
involves mainly gluons of colliding nuclei. Since the parton structure functions || point 
to a strong increase of the gluon density at small values of the below explained variable 
x, one can expect a corresponding increase of the relative contribution from charm and 
bottom decays into the intermediate mass region of dileptons with increasing beam en- 
ergy or y/s. Therefore it is a priori not obvious in which energy region one meets best 
conditions for observing a thermal dilepton signal from hot deconfmed matter. 

In this note we present a schematic view on the beam energy dependence of the 
mentioned sources of dileptons. Additionally we study systematically the influence of 
single-electron transverse momentum cuts on the invariant mass spectra of dileptons and 
show that at ALICE at LHC such cuts can suppress the Drell-Yan and correlated semilep- 
tonic decay dielectrons. Thus an observation of the thermal signal seems to be feasible, 
supposed the uncorrelated background can be accurately enough removed by like-sign 
subtraction. 



2 Beam energy dependence of dilepton sources 



Our basic equations for calculating the above mentioned dilepton sources can be found 
in Refs. [|], [J and are presented below in modified form. In Fig. 1 we show the beam 
energy dependence of the Drell-Yan yield and the yield from correlated DD, BB decays 
in central Au + Au collisions. One infers from Fig. 1 that the ratio of Drell-Yan to open 
charm (bottom) decay dileptons with invariant mass M = 2.5 GeV drops down with 
increasing yfs. Only in the region yfs < 20 (120) A-GeV the Drell-Yan yield dominates 
above the charm (bottom) decay contribution. Such a conclusion is in agreement with 
a recent analysis |J which points out that the open charm source cannot explain the 
experimental data of CERES and NA38 /NA50 of the dilepton spectra in S - U and Pb - 
Pb collisions at SPS energies. 

Also displayed in Fig. 1 are rough estimates of the thermal signal from purely de- 
confined matter assuming boost-invariant expansion. The initial parameters for RHIC 
and LHC energies are described below. For SPS energies we take the initial temperature 
Tj = 200 MeV at initial time r, = 1 fm/c, which are compatible with the observed final 
pion rapidity densities. Note that the thermal yield estimates depend sensitively on the 
assumed initial parameters. In agreement with recent findings the thermal dilepton signal 
is up to two orders of magnitude below the correlated open charm decay dileptons in a 
wide range of invariant masses at LHC, RHIC and SPS energies. Therefore, there is 
nowhere a preferred energy region; at high beam energies, however, one expects clearer 
deconfinement effects due to temperatures far above the confinement temperature. 

The results displayed in Fig. 1 cover the full phase space. Any detector acceptance will 
suppress the various sources differently [§]. Energy losses of heavy quarks propagating 
through deconfined matter |9| reduce also the open charm and bottom decay yields above 



M = 2 GeV pi 10, 11]. Notice that with increasing invariant mass the thermal yield 



drops exponentially while, for instance, the Drell-Yan yields drops less, i.e. like a power 
law. 



3 Transverse momentum cuts 

From all of these considerations the urgent problem arises whether one can find such 
kinematical gates which enable one to discriminate the thermal signal against the large 
decay background. Since the kinematics of heavy meson production and decay differs 
from that of thermal dileptons, one can expect that special kinematical restrictions su- 
perimposed on the detector acceptance will be useful for finding the needed window for 



observing thermal dileptons in the intermediate mass continuum region. In particular, 
experimental cuts on the rapidity gap between the leptons can reduce considerably the 
charm decay background |7]]. As demonstrated recently ||, the measurement of dou- 
ble differential dilepton spectra as a function of the transverse pair momentum Q± and 
transverse mass M±_ = J M 2 + Q\ within a narrow interval of M± also offers the chance 
to observe thermal dileptons at LHC In the present note we show that a large enough 
\ow-p± cut on single electrons suppresses the mentioned background processes and opens 
a window for the thermal signal in the invariant mass distribution. We take into account 
the acceptance of the ALICE detector at LHC: the single electron pseudo-rapidity is lim- 
ited to \r) e \ < 0.9 and an overall low-p_i_ cut of 1 GeV is applied (we do not impose the 
earlier planned restriction p± < 2.5 GeV). We are going to study systematically the effect 
of enlarging the \ow-p^ cut on single electrons in the invariant-mass dilepton spectra at 
yfi = 5500 GeV. 

Since the energy of individual decay electrons (positrons) has a maximum of about 
0.88 (2.2) GeV in the rest frame of the decaying D (B) meson, one can expect to get 
a strong suppression of correlated decay lepton pairs by choosing a high enough low- 
momentum cut p™ m on the individual leptons in the mid-rapidity region. For thermal 
leptons stemming from deconffned matter there is no such upper energy limit and for high 
temperature the thermal yield will not suffer such a drastically suppression by the p™ in 
cut as the decay background. To quantify this effect we perform Monte Carlo calculations 
generating dileptons from the above mentioned sources. We first employ the transparent 
lowest-order calculations and then present more complete PYTHIA simulations. 



3.1 Dilepton sources 

3.1.1 Heavy quark production and decay 

We employ here the leading order QCD processes gg — > QQ and qq — > QQ for heavy quark 
production and simulate higher order corrections by an appropriate constant K factor. 
The number of QQ, produced initially with momenta p±± = —p± 2 = P± at rapidities yi >2 
in central AA collisions can be calculated by 

dN QQ = T AA (0)JC Q H(y 1 ,y 2 ,p ± )dp 2 L dy 1 dy 2 , (1) 

f * . da Q 

H(yi,y 2 ,p±) = x 1 x 2 \fg(x u Q 2 )f g (x 2 ,Q 2 ^ 



QQ 



dt 

+ E [/ g (si, Q 2 ) f 9 (? 2 , Q 2 ) + f q (x 2 , Q 2 ) f q (x u Q 2 )] (2) 



where a®/dt are elementary cross sections (see for details 0,0]), f g , q ,q(x, Q 2 ) denote the 



parton structure functions, Xi j2 = m± (exp{±?/i} + exp{±y 2 }) / \/s and m± = yP 2 ± + Trig. 
As heavy quark masses we take m c — 1.5 GeV and m& = 4.5 GeV. We employ the HERA 
supported structure function set MRS D-' from the PDFLIB at CERN. Nuclear shad- 
owing effects are not needed to be included for our present order of magnitude estimates. 
The overlap function for central collisions is Taa(0) = A 2 /(nR A ) with Ra = 1.2A 1 / 3 fm 
and A = 200. From a comparison with results of Ref. [0] we find the scale Q 2 = ^itlq and 
JCq = 2 as most appropriate. We employ a delta function fragmentation scheme for the 
heavy quark conversion into D and B mesons. Dilepton spectra from correlated decays, 
i.e., D(B)D(B) — > e + Xe~X, are obtained from a Monte Carlo code which utilizes the 
inclusive primary electron energy distribution as delivered by JETSET 7.4. The heavy 
mesons are randomly decayed in their rest system and the resulting electrons then boosted 
appropriately. The average branching ratio of D(B) —>■ e + X is taken as 12 (10)% [Jl2"f . 

3.1.2 Drell-Yan process 

Similar to heavy quarks, our calculations of the Drell-Yan pairs are based on the lowest 
order qq e + e~ process with appropriate K factor. For the Drell-Yan production process 
of leptons at rapidities y + and ?/_ and transverse momenta = —p±- = p± we have || 

dN if Y = T AA {Q)fC DY L{y + ,y^,p ± )dp\dy + dy„, (3) 

* * - - da" 

L(y+,y-,p±) = Ys x i x 2[f<i( x iiQ 2 ) fq( x 2,Q 2 ) + f q (%2,Q 2 ) fq(xuQ 2 )\ -4-, (4) 

qq QZ 

with 

do*/dt = ^fch(y + - y_)/ (2pl[l + ch(y + - y_)] 3 ), x x , 2 = p ± (exp{±?/ + } + exp{±y_}) 
/y/s, a = 1/137, Q 2 = x x x 2 s = M 2 and K DY = 1.1. 



3.1.3 Thermal yield 

In case of thermal emission of dileptons from deconfined matter the needed distribution 
of leptons with respect to rapidities y± and transverse momenta p±± reads |8| 

dN = ^F q J drr K (^) X 2 q d 2 p ±+ d 2 p ± _dy + dy_ (5) 

with F q = Y^ q e o = | f° r u,d,s quarks, K n as modified Bessel function of nth order, and 
with dilepton transverse mass M]_ = p 2 ±+ + p\_ + 2p± + p±_ch(y + — y~). The integration is 
to be performed on the proper time r of the longitudinally expanding deconfined matter 



with temperature T(r) and quark fugacity A(r) H. Our choice of initial conditions for 



produced deconfmed matter is based on the estimates of Refs. ||, [l3j for the mini-jet 



plasma which are similar to the self-screened parton cascade model [I3|. We take as main 
set of parameters the initial temperature Tj = 1000 MeV, gluon fugacity Af = 0.5, and 



light quark fugacity Af = \Xf of the parton plasma formed at LHC at initial time = 0.2 



fm/c. For the sake of defmiteness we assume full saturation at confinement temperature 



T c = 170 MeV and a quadratic time dependence of A 9,9 (r) according to the studies [|, [15 



3.2 Results of lowest-order calculations 

The results of our lowest-order calculations of the invariant mass spectrum for various 
values of p™ m are displayed in Figs. 2-4. Comparing Figs. 2 and 3, one observes that 
already the cut p™ in = 2 GeV causes a strong suppression of the correlated charm decay 
and Drell-Yan background in the region M < 2p™ m . The above value of the invariant 
mass threshold can be estimated by using the relation M 2 = 2pj_ + p_|__[ch(y + — y_) — 
cos(0 + — </»_)], where <fr± denote the azimuthal angles of the leptons in the transverse 
plane. In order to exceed the cut p™ m most easily, the decay leptons should go parallel 
to the parent heavy mesons, which in turn are back-to-back (in the transverse plane) in 
lowest order processes. As a consequence, cos(0i — fa) ~ —1 and the minimum invariant 
mass becomes M mm pa 2p™ m for such decay pairs. For correlated bottom decay the 
electron energy is larger in the meson rest system and both leptons can easily overcome the 
threshold p™ m without such strong back-to-back correlation. Selecting, however, electrons 
with p± > p™ m = 3 GeV one can also get the corresponding threshold like behavior for 
the lepton pairs from correlated bottom decay, see Fig. 4. Therefore the thermal signal 
becomes clearly visible for such a value of p™ n due to the strong suppression of the 
considered background channels. 

The threshold behavior does not change if we include in our calculations energy loss 
effects of heavy quarks in deconfined matter || . Such effects cause mainly a suppression 
of the decay contributions (cf. 0, [10], |IT|). We also mention that shadowing effects, 
not included in our Drell-Yan and heavy quark production estimates, will diminish these 
yields somewhat. (In the mini-jet estimates of our initial conditions Q shadowing effects 
are already included). 



3.3 Results from PYTHIA 



The lowest order Drell-Yan yield has anyway M mm = 2p™ m . From the above given re- 
lations for M and M± one can derive the inequality M > 2pf n ^ 1 - (^s) • This 
relation tells us that in the region M < 5 (4) GeV only pairs with total transverse mo- 
mentum Q± > 3.25 (4.5) can contribute if p™ m = 3 GeV. Since the next-to-leading order 
Drell-Yan distribution dN / dM 2 dQ\dY drops from « to 3 - 4 GeV by nearly three 
orders of magnitude |U| one can estimate a small higher order Drell-Yan contribution in 
the small-M region. To quantify the smearing of the threshold effect by an intrinsic p± 
distribution of initial partons we perform simulations with the event generator PYTHIA 
(version 6.104 |I7||) with default switches. Results are diplayed in Fig. 5 and show that the 
sharp threshold effect from the above lowest-order Drell-Yan process is indeed somewhat 
smeared out, however, the small-M region is still clean. It turns out that the initial state 
radiation of partons before suffering a hard collision is the main reason for rising the pair 
Q± and for smearing out the sharp threshold effect, while the intrinsic p± distribution of 
partons causes minor effects. 

Let us now consider heavy quark pairs. Here, the intrinsic p± distribution, and both 
initial and final state radiations of the partons can cause a finite Q± and thus destroy the 
strong back-to-back correlation, i.e. p±i ^ p±2. In Fig. 6 we show results of simulations 
with PYTHIA for the bottom channel. The resulting primary dileptons from all correlated 
open bottom mesons are displayed. One observes that the initial state radiation causes 
a pronounced smearing of the threshold effect discussed above. Without the initial state 
radiation, as implemented in PYTHIA, the threshold effect is recovered. The intrinsic p± 
distribution and final state radiation are negligible. The conclusion of such studies is that 
an enlarged value of p™ m is necessary to keep clean the low-M region from open bottom 
decay products. For charm the smearing effect due to initial state radiation in PYTHIA 
is efficiently suppressed by the large enough low-p_i_ cut of pf m = 3 GeV. The different 
behavior of charm and bottom stem from the fact that the bottom-p^ distribution is much 
wider. As a consequence, the bottom is evolved in the average to much larger values of 
Q 2 thus experiencing stronger kicks by initial state radiation. In agreement with our 
previous findings ||, bottom therefore causes the most severe background processes at 
LHC energies. 

With PYTHIA also the dileptons from single decay chain of open bottom are acces- 
sible. This channel provides a contribution peaking at 1.5 GeV; the cut p™ in = 3 GeV 
pushes all invariant masses below 3 GeV. Therefore, unless enlarging p™ m considerably, 
it will be difficult to suppress kinematically the background below the thermal signal at 



invariant masses M < 3 GeV. Probably explicit identification and subtraction of the bot- 
tom contribution is needed in this region. Otherwise one should stress that the thermal 
signal of deconfined matter is expected to exceed the thermal hadron signal at M > 2 



GeV ||15|| . Hence the region of small values of M is not interesting in this respect. 

We conclude this subsection by mentioning that further differences between our 
above lowest-order calculations and the PYTHIA simulations stem from slightly differ- 
ent charm/bottom masses, different fragmentation schemes and decay chains. We have 
checked that for the same masses, the same fragmentation scheme (i.e. independent frag- 
mentation) and PYTHIA without intrinsic p± distribution and without initial/final state 
radiation, PYTHIA then reproduces our calculations. According to our experience, the 
stringent kinematical cuts amplify small differences in various code versions. We do not 
intent here to present a detailed prediction of the dilepton yield, but rather to demon- 
strate that an observation of a thermal signal from deconfined matter is not excluded. 
Therefore, we do not attempt any fine tuning of the codes to reproduce the bulk of data 
in pp and pA reactions. 

3.4 Physical information encoded in the continuum spectrum 

The thermal dilepton signal with single-electron low- momentum cut-offs p™ m = 2-3 GeV 
exhibits an approximate plateau in the invariant mass region 2 GeV < M < 2p™ n (see 
Fig. 4). The physical information 

encoded in the height of the plateau can be estimated by 

dN n o?R\ _ , lM , 2n /8 A „ . . x-2pf n /T l 



dM 2 dY 



4tt 2 J^/Tt \x 2 J ^J X 2 _ (M/Ti) 2 



(cf. Appendix for details of the derivation). By explicit calculation of this equation and 
as seen in Fig. 4, the height of the plateau depends sensitively on the initial temperature 
Ti and quark fugacity Af. This dependence can be used to get physical information on 
the early and hot stage of the parton matter. In particular, a variation of the initial 
temperature within the interval 0.8 - 1.2 GeV, expected as possible range of parton 
matter formed at LHC energies, leads to a considerable change of the plateau height up 
to an order of magnitude. Therefore, once an identification and a measurement of the 
thermal dilepton spectrum in the intermediate mass continuum spectrum is possible, it 
delivers some implicit information on the initial and very hot stages of the parton matter. 



4 Summary 



In summary we analyze the beam energy dependence of various expected sources of dilep- 
tons in ultrarelativistic heavy-ion collisions. Already at y/s > 20 GeV a copious produc- 
tion of charm gives rise to a dominant contribution to the dilepton spectrum at inter- 
mediate invariant mass. Taking into account the ALICE detector acceptance we study 
also systematically the effect of single-electron transverse momentum cuts. We find a 
threshold like behavior of the invariant mass spectra of dileptons from primary correlated 
charm and bottom decays and Drell-Yan yield as well: these sources are suppressed at 
M < M min ~ 2pf n for pf n > 3 GeV. In contrast to this, the thermal dilepton signal 
exhibits a plateau in this region which offers the opportunity to identify them and to gain 
information on the initial stages of deconfined matter at LHC energies. The same mecha- 
nism also works at RHIC energies, however the expected count rates are too small to make 
such a strategy feasible. The complex decay chains of heavy mesons, in particular open 
bottom, and the resulting combinatorial background make an explicit identification of the 
"hadronic cocktail" very desirable to allow a safe identification of the thermal signal. 
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Appendix 

We rewrite the general expression Eq. (|5|) for the thermal dilepton spectrum as 



= ^r F <? J drrd p+ d p- d y+ d y~ K o [-^ 



x 5 



M 2 - 2p_L + p_L_{cosh(?/ + - y_) - cos(0 + - </>_)} (7) 



and arrive after three integrations at 

dN a 2 R\ .. f , / .U, \ ,,, I 



rfM2 4?r2 F ff / dFrfMi drr^o (^J Aj(r) -I(M,M ±)P r) (8) 

with 

KM Mi - 2 1 * 2 F(<M) + A*1[E(6*)-F(^)] 

where E and F stands for the incomplete elliptic integrals 

sin</> sin0 „ „ 

,x f dt /" (1 - k z t z ) dt , , 

/ / ^(1-^(1-^2) 

with = arcsin (^^) an d A; = j^-, and A (a;) = min(x— 2p™ m , QjJ, = max(2p™ m , 
X2 = Mj_. In the special case p™ in = the expression I(M, M±,pf in ) becomes tt. In the 
small-M region one can approximate A(x) ~ Mi — 2p™ in being independent of x. Then 
the x integration yields 



l(M,M ± ,pf n ) « 2(E(7r/2,v / T^P)F((/),A;) + F(7r/2,v / r^F)E(0,A;) 

-F( 7 r/2, v / r^P)F(0,A;)). (11) 



Since for small M also Qj_/Mj_ « 1 holds, one can use E(0, 0) = F(0, 0) = and 
E(0, 1) = sin0 and finds 

Mi -2nT in 



For cuts larger than 2 GeV one gets from Eqs. (§,[T2|) the approximate expression 
Eq. (|) when approximating the r integration appropriately ||. 
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Figure 1: The dependence of dileptons from the lowest-order Drell-Yan process, and 
correlated open charm and bottom decays and the thermal source (only purely deconfined 
matter; initial conditions are described in text) on s 1//2 . Note that at SPS energies the 
hadron and a possible mixed phase (not included here) are by far the strongest thermal 
sources. 




Figure 2: The invariant mass spectra of dileptons from the Drell-Yan process, and charm 
and bottom decays, and thermal emission. The single-electron low transverse momentum 
cut is pf n = 1 GeV. 
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Figure 3: The same as in Fig. 2, but p™ in = 2 GeV. 




Figure 4: The same as in Fig. 2, but p™ m = 3 GeV and various initial temperatures (in 
units of GeV) as indicated by the labels. The fat dots indicate the estimates of the low-M 
thermal plateau according to Eq. (|]). 
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Figure 5: The invariant mass spectra of dileptons from the Drell-Yan process for p™ in = 
1, 2 and 3 GeV (from left to right). The curves depict results of PYTHIA with default 
switches. 
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Figure 6: The invariant mass spectra of dileptons from correlated open bottom meson 
decays for p™ m = 2 and 3 GeV. The curves depict results of PYTHIA (dashed curves: 
default switches, solid curves: without initial state radiation). 



